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ABSTRACT 

We propose an extension of the GC3 streamer finding method of lJohnston et all fl996) 
that can be applied to the future Gaia database. The original method looks for stream- 
ers along great circles in the sky, our extension adds the kinematical restriction that 
velocity vectors should also be constrained to lie along these great circles, as seen by 
a Galactocentric observer. We show how to use these combined criteria starting from 
heliocentric observables. We test it by using the mock Gaia catalogue of Brown et al. 
(|2005f ). which includes a realistic Galactic background and observational errors, but 
with the addition of detailed star formation histories for the simulated satellites. We 
investigate its success rate as a function of initial satellite luminosity, star formation 
history and orbit. We find that the inclusion of the kinematical restriction vastly en- 
hances the contrast between a streamer and the background, even in the presence of 
observational errors, provided we use only data with good astrometric quality (frac- 
tional errors of 30 per cent or better). The global nature of the method diminishes the 
erasing effect of phase mixing and permits the recovery of merger events of reasonable 
dynamical age. Satellites with a star formation history different to that of the Galactic 
background are also better isolated. We find that satellites in the range of 10 s — 10 9 
L Q can be recovered even for events as old as ~ 10 Gyr. Even satellites with 4 — 5 x 10 7 
L Q can be recovered for certain combinations of dynamical ages and orbits. 



1 INTRODUCTION 

The A cold dark matter (ACDM) hierarchical paradigm 
is the best model we h ave to explain the large-scale 
structure of th e Uni ve rse (ISpringel. Frenk fc White] 20061; 



missi n g satellite s prob lem (|Klypin et al .11 19991 ; I Moore et al.l 
1 19991; iKravtsov 2010 ] '). the angular momentum problem 
llAbadi et al. 1 120031; IGovernato et al. 1 12004 lOkamoto et al.l 
120051 ; IGovernato et al.l 120071 ') . or the Downsiz ing-Specific 



Spergel et aU 120071 ; iKlypin. Truiillo-Gomez fc Primackl 
Of). In recent years, it has proved to successfully repro- 
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duce a number of observational measurements (e.g. spatial 
and colour-magnitude distributions) for galaxies observed 
in the local Universe and at higher redsh ift (for recent 
reviews see iBaughl |200d ; lAvila -Reesd 120061 ) . In spite of 
the encouraging progress on the large scale, at galactic and 
subgalactic scales, the success of the model has not been 
convincingly demonstrated as yet, and a number of issues 
remain subject of a lively debate in the astronomi cal com- 
munity. Examples of these are the core-cusp issue Jde Blok 
2010] ; IGovernato et alJIgOld; iPuglielli, Widrow fc Courte~ 



Star For mation Rate relation for dwarf galaxies (|Colin et all 
20091; |Fi7mani Avila-Reese fc Rodriguez-Puebial 120091 ; 



Bauer et al. 200 



35). 



20ld ; IValenzuela et al.l 120071 ; iGentile et al.l [2004), the 
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Another powerful test for the ACDM scenario at sub- 
galactic scales, may be provided by galaxy halos, in partic- 
ular by the Milky Way (MW) stellar halo. If the Galaxy 
was assembled in a hierarchical way, there must be fos- 
sil signatures in the phase space distribution of halo stars 
and also in the Galactic disc, although kinematic features 
in the latter may have multiple origins jAntoia et al.l 120091 ; 
iMinchev et aU 120091 : iGomez fc Helmill2010l ). Theoretically, 
the origin and structure of the stellar halo have been studied 
by several authors using a variety of techniques (for a review, 
see Belmi 2008). These include cosmological numerical sim- 
ulations with and without baryonic physics |Zolotov et al.l 
2009 ; iDiemand. Madau fc Moorel 120051 ; \A~ 



lAbad i et al. 



2003) 
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and phenomenological modelling of the evolution of baryons 
inside haloes, usually in combination with N-body simu- 
lations that provide the dynamical history o f the system 
|Bullock fc Johnstonll2005l ; ICooper et alj|2009h . Historically, 
chemical and kinematic information was used as a basis to 
formulate the first galaxy formation mode l s. In t heir classi- 
cal paper, lEggen. Lvnden-Bell fc Sandagei l| 19621 ) suggested 
a rapid radial collapse that later continued to for m the stellar 
disc. About a decade later, ISearle fc Zirml (|l978l ) formulated 
the hypothesis that the stellar halo formed over a longer 
time-scale, through the agglomeration of many subgalac- 
tic 'fragments' that may be similar to the surviving dwarf 
spheroidal galaxies observed to day as satellites of the Milky 
Way. The lSearle fc Zinnl (|l978h scenario and several obser- 
vational results are in qualitati ve agreement with expec - 
tations from the ACDM model (|De Lucia fc Hemiill2 008). 
Howe ver, the observed abundance pattern ( Tolstoy et al.l 
2009) seems to exclude the possibility that a significant con- 
tribution to the stellar halo com es from disrupted satellites 
similar to the present-day dSphs ( Bullock fc Johnstonll2005l ; 
iRobertson et al.ll2005l ; front et alj|2006h . 

Over the next decade, a number of astrometric and 
spectroscopic surveys will provide accurate spatial, kine- 
matic and chemical informat ion for a large numb er of 
stars e.g. the Gaia satellite (iLindegren et al.| 20081); the 
Radial Velocity Experiment-RAVE l|Steinmetz et"aiTl2006l ) 
and the Sloan Extensi on for Galactic U nderstanding and 
Exploration- SEGUE (|Yannv et al.l l2009h . This vast fossil 
records will provide important advances in our understand- 
ing of the se quence of events which led to the formation of 
our Galaxy l|Freeman fc Bland-Hawthornll2002l ). 

It is critical to design efficient strategies in order to ex- 
tract valuable dynamical information from the plethora of 
upcoming observations. The task is not a trivial one. It is not 
clear which is the optimal strategy, or in which space sub- 
structure can be best located. Natural diagnostics like inte- 
grals of motion va riables may be hampered by the finite ac- 
curac y of surveys ijHelmi fc Whitdl 19991 : iHelmi fc de Zeeuwl 
2000). Recently, action-an gle variables have show n to pro- 
vide a promising aven ue (| Gomez fc Helmil 120101 ). In addi- 
tion, iFont et all (|2006l ) showed the usefulness of chemical 
information in identifying the relics of merger events. 

In the present work we c ontinu e the study initiated in 
iBrown. Velazquez fc Aguilarl l|2005l . hereafter B05), to re- 
visit the possibilities of extracting reliable information about 
past merger events in the halo of our Galaxy in the upcoming 
Gaia database. In the previous study, we developed the nu- 
merical machinery to build 'mock' Gaia catalogues that in- 
clude merger events from A-body simulations against a real- 
istic smooth Galactic background. Spatial, velocity and pho- 
tometric information from a Galaxy model with three sepa- 
rate components (bulge, disc and halo) was used. Sampling 
issues were tackled, like the variation in probing depth of the 
satellite stellar luminosity function, as a function of position 
along a tidal streamer. A realistic model for Gaia errors was 
included too. With the resulting c atalogue, a re-exam i nation 
of the integral-space method of IHelmi fc de Zeeuwl (|2000l ) 
was done and it was concluded that, although promising, 
in real practice observational errors and the presence of an 
overwhelming Galatic background can compromise the ef- 
ficiency of this method and great care in its use must be 
exercised. It is clear that we must use as many search tech- 



niques as possible, preferably those that include as much 
observationally available information as possible. 

Following this path, we have decided to investigate now 
the 'Great Circle Cell Count ' (GC 3) method proposed by 
Ijohnston. Hernquist fc Boltd l|l996l) . in order to find a way 
to increase its efficiency. We propose an extension to this 
method, that we dubbed mGC3 ('modified GC3'), which 
adds kinematical information to the original method, in or- 
der to detect substructures in the MW stellar halo. We have 
used the mock Gaia catalogue of B05, with the further im- 
provement of a more detailed modeling of the simulated pho- 
tometry of the satellites, so they now reflect specific histories 
of star formation. The expected Gaia measurement errors 
were simulated as a function of magnitude, colour, and sky 
position of the stars, as described in appendix A of B05. 
The values of the expected errors were updated to the latest 
predictions of the scientific performance of Gaia. The lat- 
ter can be found at http://www.rssd.esa.int/gaia under the 
'Science performance' entry. Compared to the performance 
predictions used in B05 the current astrometric error esti- 
mates are larger by a factor of two. We specifically asses 
the effect of the predicted Gaia observational errors on the 
ability of GC3/mGC3 strategies to find simulated satellite 
remnants. We also study the effects on recovery efficiency 
when varying the orbit, luminosity and star formation his- 
tory of the satellites. 

In section[2l we present a brief outlook on the problem of 
identifying past merger events and review the original GC3 
method. Some specific examples are shown to illustrate how 
it works. In section[3]the mGC3 extension is introduced. The 
basic concept, as well as a practical way to implement it, are 
discussed. In section [4] we evaluate the applicability of the 
mGC3 method using the mock catalogue of B05. We also 
describe the way in which we have included various star for- 
mation histories for the satellites using a nd adaptation of the 
stellar population synthesis software of iBruzual fc Charlotl 
(2003). In section [S] we study in detail the efficiency of the 
mCG3 method in identifying past merger events, as a func- 
tion of satellite orbit, luminosity and star formation history. 
Our conclusions are presented in section [5] 



2 THE GREAT CIRCLE CELL COUNT 
METHOD 

Current computer capabilities allow us to explore the forma- 
tion and evolution of stellar streams under a variety of condi- 
tions, such as different orbits and underlying potentials, via 
N-body simulations. The understanding of common proper- 
ties of stellar streams, as well as their time evolution, are key 
to develop effective identification techniques. This goes hand 
in hand with the observational developments that will allow 
the practical application of such techniques to real data. 

Key properties that determine the dynamical structure 
and evolution of stellar streams, which in turn can be used 
in their identification, are the integrals of motion: the total 
energy E, if the Galatic potential is stationary; the total an- 
gular momentum L, if the potential has spherical symmetry; 
and its projection L z , if the potential has axial symmetry. 

The existence of these conserved quantities leads natu- 
rally to the search of strea ms as entities in int e grals of mo- 
tion space, as proposed by IHelmi fc de Zeeuwl fcoool ). This 
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Figure 1. Illustration of spatial and corresponding GC3 pole count maps of simple stellar distributions, shown in Aitoff projections, (a) 
Spatial distribution of 9 X 10 3 simulated stars in a planar circular orbit, with an inclination of 45° with respect to the Galactic Plane. 
(b) GC3 pole count map of stars shown in (a), (c) Spatial distribution of 9 X 10 3 simulated stars with a localized gaussian distribution. 
(d) GC3 pole count map of stars shown in (d). 



technique requires the computation of the conserved quati- 
ties in terms of observed ones, like (I, b, w, v r , Hi, Hb), the 
Galactic latitude and longitude, parallax, radial velocity and 
proper motions, respectively. These computations involve 
non-linear transformations resulting also in non-linear error 
propagation which, as B05 show, lead to systematic errors 
in E, L, L z as their expectation values do not equal their 
true values once the observational errors are introduced. 
Furthermore, the spread of errors is greatly magnified, as 
shown there. Additionally, the computation of the total en- 
ergy E requires an assumption about the underlying po- 
tential, which is an unknown and thus introduces another 
source of error. It is clear that, although powerful in theory, 
the integrals of motion technique has its limitations when 
applied to real data. It is thus convenient to look for addi- 
tional search techniques that may complement it. 

A different approach in the identifi cation of tidal 
stream s is the GC3 method proposed by Ijohnston et al.l 
( 1996), which uses the conservation of the total angular mo- 
mentum in a geometrical way. It relies on the fact that the 
orbits of stars in a stream in a spherical potential will be 
confined to a plane that contains the potential center. As 
seen from the Galactic center, this plane is a great circle on 
the sky. In the more general case of an axisymmetric Galac- 
tic potential, the orbital plane precesses around the axis 
of symmetry of the system. This will result in the stream 
being approximately confined to a plane with a finite thick- 
ness, which in turn corresponds to a great circle of finite 
width. A further complication arises from the fact that we 
are not at the Galactic center, which spoils to some degree 
the great circle effect, particularly for streamers whose orbits 
lie within, or near the Sun's orbit. 



The practical implementation of the GC3 method be- 
gins by laying down a grid on the celestial sphere. Each grid 
point is then associated to a unit vector L (pole vector) that 
points from the observer to the grid point. All stars on the 
sky that lie on a band defined by the great circle in the 
plane orthogonal to the pole vector, are associated to this 
pole or grid point, the angular width of the band being spec- 
ified independently. The number of stars within this band 
are assigned to the corresponding pole. The counts associ- 
ated to every grid point in the mesh thus define a pole count 
map. Significant maxima in this pole count map are then ac- 
cepted as possible streamers. Because of the geometry, the 
pole vector will be parallel to the angular momentum of all 
stars in the band whose velocity vector also lies in the or- 
bital plane defined by the great circle. This latter property 
will be exploited later on, when we introduce the extended 
GC3 method. 

|L ■ f| sC 5 r (1) 

where f is a unit vector pointing to the star and S r — 
sin 56 is the tolerance which allows for the width 80 of each 
great circle associate d to a cell (see Sec. 3.1 and Fig. 7 in 
Ijohnston et al.l 1996 V The distribution of pole counts is then 
represented in a map of the celestial sphere. 

We illustrate the pole count maps for two simple sim- 
ulated stellar distributions in Fig. []J,-d, in order to explain 
their morphology. The spatial distribution of 9 x 10 3 stars 
in a circular orbit with a 45° inclination is shown in Fig. [1J, 
and the corresponding pole count map in Fig.[]J>. This shows 
how the signature in a pole count map, of stars in a planar 
orbit, corresponds to a localized peak which coincides with 
the direction of the normal vector defining the orbital plane. 
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Since two antiparallel pole vectors define the same plane or 
great circle in the celestial sphere, the information on one 
hemisphere of a pole count map is replicated on the other; 
hence the two peaks at exactly opposite directions in Fig[lp. 
Figure []J shows a map of stars in a localized stellar distri- 
bution, with the corresponding pole count map in Fig. [TJi. 
The maxima in pole counts are in this case distributed along 
a great circle. This happens because all the possible planes 
which contain both, the localized distribution of stars and 
the observer at the Galactic Center, have poles which are in 
turn contained in a specific plane, of which the projection 
on the sky is the great circle seen in Fig. [TJi. 

In the GC3 method it is expected that the stars in the 
Galactic background will contribute a smooth distribution 
across the entire map, modulated only by the stellar density 
dependence on position in the sky. The stars in a stream, 
on the other hand, will contribute only to the cells which 
coincide best with the stream's orbital plane, creating a local 
maximum in the pole count distribution. 

The GC3 me thod was used on a C-star survey by 
llbata et alj (|200ll ) in order to quantify the statistical sig- 
nificance of a great circle C-star overdensity and to identify 
it as part of the tidal stream from the Sagittarius (Sgr) 
dSph galaxy. Lat er, it was used on 2MASS M-giants by 
llbata et al.l (|2002l ), on a search for tidal debris in the Galac- 
tic Halo. They succesfully recovered the Sgr stream but 
failed to indentify any other significant features. 

In the following subsections we explore the performance 
of the original GC3 method by applying it to a simulated, or 
mock Gaia catalogue (B05). The catalogue includes obser- 
vational uncertainties and a realistic number of stars in the 
Galactic background. It illustrates the need to modify the 
original method in order to reduce the contribution of the 
Galactic background, and increase its sensitivity to detect 
lower luminosity tidal streams. 



2.1 The Mock Gaia Catalogue 

The Galactic catalogue from B05 constitutes a random re- 
alization of kinematic properties of Milky Way stars in the 
Galactic disc, bulge and halo, with their corresponding den- 
sity, velocity, age and color distributions, as well as the ap- 
propriate normalization in total luminosity. This resulted in 
a mock Gaia catalogue containing ~ 3.5 x 10 s observable 
Milky Way stars with full phase-space information and ob- 
servational errors as are expected from Gaia, as well as the 
error-free quantities for comparison. The two most impor- 
tant aspects to be emphasized about the mock catalogue are 
the realistic simulation of both the observational errors and 
the total number of stars that will be observable with Gaia, 
which will be key issues in determining the applicability of 
any stream searching method. 

2.2 GC3 Galactic pole count maps 

We restricted the catalogue to stars with |6| > 10°, in order 
to avoid the Galactic Plane. The tolerance used in {JJ was 
59 — 5° which corresponds to the half-width of each great 
circle cell and the GC3 pole counts were computed on a 
72 x 72 cell grid, uniformly spaced on the surface of the 
celestial sphere. We chose the tolerance to be slightly less 
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Figure 2. GC3 pole count map of Galactic background from 
error-free mock Gaia catalogue, in a Galactoccntric reference 
frame. The color scale represents the number of stars per pole 
as shown in the color bar. 



than the 50 = 6° half-width which llbata et all (|2002h find 
to maximize the signal to noise ratio of the Sgr tidal stream 
feature in their 2MASS M-giant pole count maps. 

We computed the pole counts for the error-free mock 
Gaia catalogue, using the GC3 method's position criterion 
expressed in {JJ in a Galactocentric reference frame, i.e. 
position vectors and pole coordinates are Galactocentric, 
the latter corresponding to normal vectors of planes that 
go through the Galactic center. The resulting pole count 
map is shown in Fig. [2] In this reference frame we name 
the longitude and latitude angles <f> and 9 respectively, with 
the Galactic Plane at 9 = 0°, the North Galactic Pole at 
9 — +90° and cf> — 0° in the direction away from the Sun. 

In this reference frame the pole count map would be 
expected to be uniform in <f) and with a positive gradient in 
9 (minimum at the equator) , because of the axial symmetry 
of the Galaxy about the Galactic center and the latitude 
dependence in star density. However, the pole count map 
shown in Fig. [5] clearly does not follow this pattern in <j>. 
This is due to both the exclusion criterion in Galactic lati- 
tude (|6| > 10°) and the fact we are considering only stars 
which will be observable by Gaia; both criteria being in- 
herently heliocentric. First, the exclusion criterion imposed 
on the heliocentric Galactic latitude does not filter stars in 
the Galactic Plane near the Sun at high latitudes, these stars 
contribute to pole counts for poles in the redish circle in Fig. 
[2] with their maximum contribution being for <j> P oie ~ 90° 
and 4> P oie ~ 270°. In addition, planes with 4> po i e = 0°, 180° 
do not contain solar neighbourghood stars and thus include 
stars at heliocentric distances that are larger on average than 
for planes that do go through the Sun with 4>poie = 90° , 270° 
(this is illustrated for 4> po ie ~ 0° and (f> vo ie ~ 90° in Fig. [3}. 
Therefore the number of observable stars that contribute to 
poles 4>poie ~ 90° , 270° will be larger than in the perpendic- 
ular direction. These two effects reinforce one another and 
give rise to the morphology of the pole count map of Fig. [2] 

The GC3 method's capability for detecting tidal 
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Figure 3. Projections of mock Gaia catalogue stars. The color 
scale indicates heliocentric distances Dh e ii as indicated by the 
color bar. Top: X-Z projection. The grey shaded area shows the 
plane with pole <j>pole = 0° and 8 = 0°, which does not con- 
tain solar neighborhood stars. Bottom: Y-Z projection. The grey 
shaded area shows the plane with pole <j> po i e = 90° and 8 = 0°. 
This plane contains solar neighborhood stars. Also, as shown by 
the color scale, the heliocentric distances of stars in the shaded 
area of the top panel are on average larger than those in the 
shaded area in the bottom panel. These effects give rise to the 
morphology noted on Fig. [2] (see text for details). 



streams clearly depends on how large is the contribution of 
the stream with respect to the Galactic background. In Fig. 
[5] the mode of pole counts is ~ 4.5 x 10 6 stars/pole, which 
is a measure of the typical contribution of Galactic stars in 
this error-free map. On the other hand, a dwarf galaxy with 
a stellar population similar to the Galactic Halo (i.e. an age 
of 13.6 Gyr and iron abundance [Fe/H] = —1.7) and total 
luminosity of 10 7 Lq has ~ 1.2 x 10 6 stars brighter than 
Mv = 5, according to the stella r population synthesis mod- 
els of iBruzual fc Chariot! (|2003h . The number of observable 
stars, however, will depend on the distance distribution of 
the stars in the stream, but will necesarily be smaller, spe- 
cially since most of these are rather faint (My > 4) main- 
sequence stars. Therefore, the signature of such a dwarf 
galaxy in a pole count map would be several times below 
the typical number of Galactic background stars and thus 
very hard to detect. 

As exhibited in this example, the number of Galactic 
contaminants in a cell is an important limitation of the GC3 
method. In the following section we propose a modification 
of the method, to help discriminate stream stars from the 
smooth halo using the main attribute of the original GC3 
method, namely the use of angular momentum conservation 
through a geometrical approach. 



3 THE MODIFIED GC3 METHOD 

The GC3 method of Johnston et al., as explained above, 
uses solely positional information. In order to improve the 
'signal' of a stellar stream over the 'noise' of the smooth 
background halo population, we propose the inclusion of ve- 
locity information as well. Since the stellar stream is roughly 
confined to a plane that contains the center of symmetry of 
the potential, as seen from this center both its position and 
velocity vectors will be contained in the orbital plane. 

A grid of all the possible great circle cells is devised, 
with each cell uniquely determined by its pole, that is the 
unit vector L which is perpendicular to the plane it defines. 
The modified GC3 (mGC3) method proposed here for count- 
ing the number of stars associated with each pole, can be 
expressed in terms of the following Galactocentric position 
and velocity criteria 



|L ■ f I < 5 r and L ■ v| ^ 5 V 



(2) 



where f and v are unit vectors parallel to the star posi- 
tion and velocity vectors, 5 T and S v are the tolerances which 
allow for the width of each great circle associated with a 
cell. The pole vector L that corresponds to the cell which 
best coincides with a stream is parallel to the streams's total 
agular momentum. 

As with the original GC3, the mGC3 criteria in ((2| 
will hold for tidal streams evolved in a spherical or slightly 
flattened potential. For the Milky Way (MW), according 
to several authors, the orbit of the Sgr tidal stream con- 
strains the inner halo (< 60kpc) to be slightly flattened 
with a Z-axis flattening q z = 0.85 — 0.95 for oblate models 
( Johnston. Law fc Maiewski 1120051 : iMartmez-Delgado et al.l 
12004 llbata et alll200ll), and q z = 1.25 for prolat e mod- 
els l|Helmill200i; lLaw. Johnston fc Maiewski 1120051 ); albeit 
more recently lLaw. Maiewski fc Johnston I (|2009l ) find that 
a triaxial halo with q z = 1.25 and q y = 1.5 provides a better 
fit of the radial velocity and distance distribution of ob- 
served Sgr stream stars. On the other hand, depending on 
dynamical age, orbital inclination and how close the stellar 
system comes to the Galactic plane, the disc's potential can 
perturb the symmetry of the overall potential felt by the 
stream. Ijohnston et al.l l|2008l ) classify stream morphologies 
resulting from their N-body simulations in three categories 
'great circle', 'cloudy' and 'mixed', illustrated in their Fig. 
1, as well as transition types in between these categories, 
as shown in their Fig. 2. The mGC3 criteria proposed here 
will hold for the 'great circle' streams, as well as for the dy- 
namically younger parts of 'cloudy' and 'cloudy-great-circle' 
streams, and clearly not for 'mixed' morphology streams. 



3.1 Practical Implementation of the mGC3 
criteria 

The position and velocity criteria proposed above, need to 
be written in terms of observable quantities. Our first step is 
to write them using the Galactocentric position and velocity 
vectors of each star: 



|L • r ga i| ^ ||r ga i||<5r and |L • v ga i| 5? ||v ga i||<5„ 



(3) 



where 8 r = sin<5^v, S v = sin<5?/>„ and Sip r ,8ip v are 
the complements of the angles between L and the vectors 
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r ga i, and v ga i respectively, which correspond to the toler- 
ance width of the great circle associated with the cell. 

The vectors r ga i and v ga i are in turn expressed in terms 
of the observable quantities (l,b,U7,v r , fix, fib) as follows 

r gai = T0 + Apzu" 1 [(cos I cos i>)x + (sin / cos b)y + (sin 6)zj 

Vgai = Vq + V r f + ZU^ 1 \(A,,tU COS 6)1 + (Aj,it{,)bJ , 

where A p = 10 3 mas-pc, A v = 4.74047 yr km s" 1 ; 
{x, y, z} are the unit vectors in the cartesian Galactocentric 
reference frame and {f, 1, b} are the unit vectors in a spher- 
ical heliocentric reference frame. The latter depend upon 
(l,b), though not on the parallax w, when expresed in the 
cartesian Galactocentric reference frame (see Appendix [X]). 

As can be seen from ([4]), the expressions for the Galac- 
tocentric position and velocity vectors r ga i and v ga i de- 
pend upon the reciprocal of the parallax, which as dis- 
cussed above, will introduce systematic errors in the deter- 
mination of both vectors. Therefore we define the vectors 
r'g a i = ccrr ga i, and v' ga i = ra7V ga i. With these definitions the 
parallax will enter multiplying instead of dividing (see [5]). 

r'gai = zutq + j4 p (cos I cos 6)x + (sin / cos b)y + (sin 6)z 
v'gai = mv Q + zuv r r + (A v (j,i cos 6)1 + (A v fi b )h 

We can thus rewrite the criteria of © in terms of the 
vectors r' ga i and v' ga i as shown in which will constitute 
our working criteria. 

|L • r'gai | < ||r' g ai||5r- and |L ■ v' ga i| < ||v' ga i||6„ (6) 

These expressions are numerically equivalent to those in 
(0), though they are in practice more useful, since the effect 
of error propagation due to the reciprocal of the parallax 
is avoided. The method based on these criteria relies on an 
accurate knowledge of and Vq, which we can expect to 
be accurately determined from Gaia data. 

3.2 Modified GC3 applied to the mock Gaia 
Catalogue 

We computed the pole counts for the mock Gaia catalogue 
described in Sec. 12.11 now using the modified GC3 (mGC3) 
criteria expressed in ([6]). As in Sec.[2l we restricted the cat- 
alogue to stars with \b\ > 10° and the mGC3 pole counts 
were computed on a 72 x 72 cell grid, with equal tolerances 
for the position and velocity criteria fo/v = 6ip v — 5°. 

3.2.1 mGC3 pole counts of mock catalogue without errors 

For the error-free mock catalogue the resulting pole count 
map is shown in Fig. 2] In this map the mode of pole counts 
is ~ 4 x 10 5 , a factor of 10 lower than the typical pole counts 
in the GC3 map of Fig. [2]This is a consequence of the veloc- 
ity requirement of which decreases the overall number 
of stars per pole because the contribution to any given pole, 
or great circle cell, does not come from all the stars in the 
plane of the great circle associated with the cell, but comes 
only from those whose velocity is also contained in the great 
circle plane, therefore decreasing the contribution to each 
pole by chance alignments. This effectively lowers the pole 
counts on the whole map and also reduces the strength of 
the pattern of excess pole counts around (j> po ie = 90° , 270° 




2e+05 3e+05 4e+05 5e+05 6e+05 7e+05 8e+05 9e+05 le+06 
Nstars (#stars/pole) 



Figure 4. Map of mGC3 pole counts for mock Gaia catalogue 
without observational errors. The color scale represents the num- 
ber of stars per pole as shown in the color bar. Notice the change 
in normalization with respect to Fig. [2] 
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Figure 5. Map of mGC3 pole counts for mock Gaia catalogue in- 
cluding simulated observational errors. The color scale represents 
the number of stars per pole as shown in the color bar. Note that 
the maximum of the color scale in this figure is ~ 10 per cent of 
Fig.H 

because, as discussed in Sec. [2j the majority of stars that 
contribute to this excess are in the Galactic Plane, and so 
their velocities do not lie on the plane associated to these 
(fipoie values. 

3.2.2 mGC3 pole counts of mock catalogue with errors 

The pole count map for the mock catalogue including the 
simulated errors is shown in Fig. [5] and was computed using 
only stars with |6| > 10° and with positive parallax with 
errors smaller than 30 per cent. 

This requirement on the parallax causes a pattern of 
excess pole counts for poles at <f) vo i e ~ 90°, 270°, similar to 
that seen in the GC3 pole count map discussed in Sec. [2] 
The origin of this pattern is essentially the same as discused 
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Table 1. Orbital parameters of simulated satellites 



Sat 


Mass 


Pericentre 


Apocentre 


Inclination 


ID 


(xlO 7 M ) 


(kpc) 


(kpc) 


angle (°) 


SI 


5.6 


8.75 


105 


30 


S2 


5.6 


7.0 


60 


15 


S3 


2.8 


7.0 


80 


60 


S4 


2.S 


40 


60 


25 


S5 


2.8 


3.5 


55 


45 



in Sec. [2] great circle cells with poles 4>poie ~ 0°, 180° include 
stars with heliocentric distances that are on average larger 
than for perpendicular cells; since parallax errors increase 
with heliocentric distance, this causes less stars to pass the 
parallax error filter for cells that do not include the Sun, and 
conversely an excess for cells that do (4>poie ~ 90°, 270°). 
However, the background level is vastly reduced compared 
to that of Fig. H 



4 APPLICABILITY OF THE MGC3 METHOD 

We now evaluate the applicability of the mGC3 method by 
applying it to pole count maps including the Galactic back- 
ground from the mock Gaia catalogue described in Sec. 12.11 
and tidal streams from N-body simulations also from B05. In 
the following section we describe the kinematic and photo- 
metric characteristics used in our experiments and the mor- 
phology of mGC3 pole count maps is illustrated in detail in 
Sec. 14.21 for a typical tidal stream. 

4.1 Assigning observable properties to simulated 
satellite streams 

We use the simulated satellites from B05 (see paper for de- 
tails), which are the results of N-body simulations of self- 
consistent systems with 10 6 particles evolved in a rigid ax- 
isymmetric Galactic potential. The Galactic potential used 
was oblate with a flattening in the potential of q$> h — 0.8. 
The satellites have two different initial masses of 2.8 x 10 7 
Mq and 5.6 x 10 7 Mq and different pericentre, apocentre and 
inclination angles with respect to the Galatic plane, which 
are summarized in Table [T] The satellite IDs used here (Sl- 
S5) correspond to the Run No. (1-5) used in Table 6 of B05. 
The resulting orbits are illustrated in Fig. 4 in B05. For 
each satellite, snapshots of the distribution of N-body parti- 
cles are available for dynamical ages ranging from to 9.75 
Gyr at 0.65 Gyr steps. We use the term dynamical age to 
emphasize this does not refer to the age of the stellar pop- 
ulation, but to the elapsed time in the N-body simulation. 
We take the dynamical age to be zero in the beginning of 
the N-body simulation, when the satellite is entirely bound. 

Embedding the simulated satellite streams in the Galac- 
tic background described in Sec. 12.11 requires that we make 
a realistic simulation of the number of stars that will be 
observable by Gaia in each of the simulated satellites. This 
number depends on the total luminosity Ly of the satellite; 
on its orbit, which determines the distribution of observer 
to star distances along the stream; and on its star formation 
history (SFH), which will determine the luminosity function 
of the stellar population at any given age. 



For our experiments we varied the total luminosity Lv 
in the range 1 x 10 7 Lq to 1 x 10 9 Lq. In all cases the 
luminosity Lv refers to the total luminosity of the stars in 
the whole satellite system including the tidal tails, not just 
the bound core. 

For each simulated satellite we assigned two different 
SFHs, whose properties are described in Table [2] These are 
termed Halo-type and Carina-type since they have been tai- 
lored to resemble the SFHs of the Galactic Halo and Ca- 
rina dSph satellite, respectively. As can be seen from Ta- 
ble [21 the Halo-type SFH consists of only an old burst 
with an age of 13 Gyr, with [Fe/H] near the peak of 
the iron abund ance distribution of halo stars ([Fe/H]=-1.8, 
|Prantzosi r2009). The Carina-type SFH consists of three sep- 
arate bursts producing old (13 Gyr), intermediate-age (8 
Gyr) and young (3 Gyr) populations, the intermediate and 
old being the dominant ones, and a single metallicity (al- 
though Carina exhibits a range in [Fe/H] from -2.0 t o -1.0, 
see e.g. ICarigi fc Hernandez! 120081 ; ICarigi et al.ll2002l ). This 
choice of SFHs and [Fe/H] does not intend to accurately 
represent the stellar populations of neither the present-day 
dSphs nor of the possible halo building-blocks. It only in- 
tends to illustrate the performance of the mGC3 method in 
the least favorable case (the Halo-type SFH) and a slightly 
more favorable one (the Carina- type SFH), in terms of the 
number of stars observable by Gaia. 

Using an adaptation of the iBruzual Sz Charlotl (|2003l ) 
stellar population synthesis software, we generated a ran- 
dom realization of the photometric properties of stars in 
a system with total luminosity Lv and a given SFH and 
age. This is done by selecting the SFH, and metallicity Z 
(a metallicity Z=0.0004 corresponds to an iron abundance 
[Fe/H]=-1.7 at a solar [a/Fe] ratio). The ages of the stellar 
populations were fixed at the values given by the SFHs in 
Table [2l independently of the dynamical age, for all exper- 
iments. Luminosities and temperatures are then randomly 
drawn for as many stars as necessary, until the total lumi- 
nosity reaches the desired value of Ly. Since a large number 
of the stars in the realization will be low-mass stars too faint 
to be observed, we keep only the photometric information of 
stars brigther than an arbitrary cutoff at absolute magnitude 
Mv = 5. Although we discard these very faint stars, it is im- 
portant to emphasize that their luminosities do contribute 
to the total luminosity of the system. The magnitude cutoff 
was chosen at Mv = 5 since it is fainter than the turn-off of 
a 13.4 Gyr old population of halo-like metallicity (My ~ 4). 
This ensures that the turn-off will be brigther than our cutoff 
for a population of any given age. Nevertheless, the fraction 
of these stars that will be observable will be determined by 
the Gaia magnitude cut-off (V < 20). Following this recipe 
ensures that our random realization accurately represents 
both the luminosity function of the stars given by the SFH, 
as well as the total number of stars in a system with total 
luminosity Lv- 

Finally, we randomly assign the photometric proper- 
ties generated with this procedure, to the N-body particles 
of the simulated satellites as explained in Sec. 4.3 in B05, 
and then compute, using the recipe detailed in Sec. 4.1 in 
B05, the Galatic longitude, latitude, parallax, radial veloci- 
ties and proper motions with their corresponding errors, for 
each simulated satellite star within the completeness limit 
of Gaia (V < 20). 
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Table 2. Star Formation Histories 



Name 


Age a 


[Fe/H] 


Lv 




(Gy~ 


(dex) 


(per cent) 


Halo-type 


13 


-1.7 


100 




13 


-1.7 


30 


Carina-type 


8 


-1.7 


50 




3 


-1.7 


20 



a Age of the star formation bursts assigned to simulated dwarf galaxie 

4.2 Aplication to a typical stream 

In this section we use a typical satellite stream with a chosen 
dynamical age, luminosity, orbit and SFH in order to explain 
the morphology of the pole count map and the methodology 
devised to detect in it the signature produced by the stream. 
For the present and all following experiments, pole count 
maps were produced with a 72 x 72 spherically uniform grid, 
position and velocity tolerances of 5ip r = 5tp v = 5° and in- 
cluding from the simulated catalogues only those stars with 
positive parallax with errors less than 30 per cent and with 
\b\ > 10°, to avoid the Galactic Plane. 

For this example we chose satellite S2 at a dynamical 
age of 5.85 Gyr, with an elongated and inclined orbit, a 
Carina- type SFH and total luminosity of 3 x 10 8 Lq. The 
corresponding Galactocentric sky distribution is shown in 
Fig. [5] The color scale of the figure indicates the particles' 
Galactocentric distances. The mGC3 pole count map which 
corresponds to this simulated stream alone is shown in Fig. 
[7^; and the mGC3 pole count map of this stream embedded 
in the Galactic background (with errors) described in Sec. 
13.2.21 is shown in Fig.[7p. In this figure, the signature of the 
stream in pole counts is a barely perceptible excess around 
WVie,6» poie ) = (85°,35°fl This is the kind of feature we 
need to detect in pole count maps in an automated fashion, 
in order to evaluate the method's capabilities to detect tidal 
streams of various characteristics. 

Since the Galactic background leaves a smooth signa- 
ture in the pole count maps, we can use the standard image 
processing technique of unsharp masking to remove its con- 
tribution. Unsharp masking consists in subtracting from the 
original image a smoothed image, in which the value of each 
pixel corresponds to the median value in its neighbourghood. 
The subtracted image or, in this case, the subtracted pole 
count map has a much more uniform background and the 
contrast of localized excesses is enhanced. The smoothed 
image, made from the pole count map on Fig. [7p using a 
neighbourhood of 5 x 5 pixels, is shown in Fig. [7J;. The sub- 
tracted pole count map is shown in Fig. (TJl for the present 
example. In this subtracted map the stream's signature is 
much more evident than it was in the original pole count 
map (Fig. [7p). The color scale in the subtracted map indi- 
cates the amplitude or height of the excesses in units of the 
background's standard deviation a, and clearly shows that 
the stream in this example is detected at the ~ 7a level in 

1 The slight differences in the shape of map features with 9 < 
0° compared to those with 8 > 0° are an effect of the Aitoff 
projection. 




40 
Rgal [kpc) 



Figure 6. Galactocentric sky distribution of N-body particles 
in satellite S2 at a dynamical age of 5.85 Gyr. The color scale 
represents Galactocentric distance as shown in the color bar. 



the subtracted pole count map; whereas in original map the 
excess appears only at the ~ 2cr level. 

The subtracted pole count map of Fig. [7jl still shows 
signs of non- uniformity. Some regions in the subtracted map, 
particularly near an excess, appear with pole counts below 
the background (i.e. negative sigma amplitudes in Fig. 01). 
This is caused by an oversubtraction of the background due 
to the excess itself, which increases the median value of the 
pole counts in pixels surrounding it in the smoothed im- 
age. Also, spurious detections tend to arise at the 2 — 3.5cr 
level due to imperfections in the background subtraction. In 
the examples of following sections we consider bona fide ex- 
cess detections as those corresponding to excess counts with 
amplitudes larger than 4er. We also require |# po ; e | < 80°, 
consistently with our avoidance zone of |b| > 10°. 

In the following section we systematically apply this 
procedure to recover the simulated streams from the detec- 
tion of excesses in mGC3 pole count maps. 



5 VARIABLES AFFECTING THE 

DETECTABILITY OF TIDAL STREAMS 

In this section we explore the efficiency of the method in 
recovering simulated streams under different conditions of 
total luminosity, SFH and orbital parameters. 



5.1 Total Luminosity and Star Formation History 

The effect of increasing the total luminosity Lv of a system 
is to increase the number of observable stars proportionally. 
The stream's signature increases with total luminosity, as 
can be seen in Fig. [H] For all five satellites, at an arbitrarily 
fixed dynamical age of 7.15 Gyr, Fig.|8]shows the amplitude 
of the excesses recovered as a function of Lv- The signal 
to noise ratio (SNR, understood as height of the peak over 
mean background, in terms of the background's a) of the 
stream increases with Lv- For some satellites, i.e. S2, the 
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Figure 7. (a) Map of mGC3 pole counts for stars in satellite S2 alone, at a dynamical age of 5.85 Gyr. (b) Map of mGC3 pole counts 
for satellite S2 (panel a) embedded in the mock Galactic background, (c) Smoothed image of the pole count map in panel (b) used for 
the 'unsharp masking'. The color scales of panels (a), (b) and (c) represent the number of stars per pole as shown in the color bar. Note 
that the range of the color scale of panel (a) is 5 times smaller than that of panels (b) and (c). (d) Subtracted pole count map. The color 
scale and contours represent the amplitude or height of the excesses in units of the background's standard deviation cr. 



SNR reaches a plateau since, depending on the shape of the 
stream, the pole count signature can differ from a localized 
peak if the stream is very disrupted (see stream orbits in 
Fig.4 of B05). Therefore increasing the luminosity can also 
increase the dispersion of pole counts in the background, 
maintaining the height approximately constant when ex- 
pressed in terms of the background's standard deviation, 
as in Fig. [8] 



Figure[8]also shows the effect of varying the SFH, which 
in terms of the pole count maps is only reflected in a change 
in the fraction of observable stars. At the same total lumi- 
nosity, a SFH having young or intermediate-age bursts (i.e. 
a Carina-like SFH) yields a larger fraction of brigther stars 
compared to a population with only an old burst (i.e. an 
Halo-like SFH). This can be seen in Fig. [8] were for most 
satellites at a given luminosity a stream is detected at a lower 
sigma level if the population is given by the Halo-like SFH 
(triangles) than if it is given by a Carina- like SFH (circles). 
For very disrupted satellites like S2 and S5, as mentioned 
before, the signature in pole count maps differs from a lo- 
calized peak; therefore the significance of the detection satu- 
rates earlier in luminosity for brighther populations causing 
the stream to be detected at higher sigma levels for fainter 
populations. 



Finally, it can be seen that satellite SI can be detected 
with the mGC3 method down to a luminosity as low as 
~4x 10 7 L.0, even for this old dynamical age of 7.15 Gyr. 
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Figure 8. Height of detected maxima on pole count maps as a 
function of the satellite's initial total luminosity Ly. 



5.2 Orbital parameters and dynamical age 

Orbital parameters are clearly important since they deter- 
mine the shape of the orbit. How disrupted the stream will 
become strongly depends on the perigalacticon and incli- 
nation angle. Nevertheless, the most influential parameter 
affecting the detectability is rather the distribution of helio- 
centric distances of the stream stars at different dynamical 
ages, since this will determine the fraction of stars that will 
be observable in our simulated survey. 

For all dynamical ages of the five satellite's streams, 
we computed -Dheiio the mean of the heliocentric distance 
distribution, as a measure of a representative heliocentric 
distance of the system. The plot of Dhciio as a function of 
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Figure 9. Heliocentric distance -Dhelio as a function of dynamical 
age for all five satellites. 



dynamical age in Fig. shows that as the dynamical age in- 
creases, stream stars can come closer to the Sun as they are 
gradually distributed along the orbit, hence decreasing the 
mean of the distribution of Dhoiio- This effect would be en- 
hanced by the action of dynamical friction which, although 
not included in the simulations, is expected to contribute in 
the evolution of real systems. 

Therefore, the detectability of streams should increase 
with decreasing Dhelio- This is illustrated in Fig. 1101 a plot 
of the a detection level of the simulated satellite streams 
in pole count maps, as a function of Dhdio, for satellites 
with a Carina-type SFH and two different luminosities Lv = 
2 x 1O 8 L (bottom) and L v = 2x 1O 9 L (top). These plots 
show that old streams (> 7 Gyr) at heliocentric distances up 
to .Dhoiio ~ 30 kpc and Dhelio ~ 40 kpc can be detected at a 
high sigma level, for the bright (Lv — 2 x 1O 9 L ) and faint 
satellites (Lv = 2 x 1O 8 L ) respectively. On the other hand, 
streams can be detected at larger distances only if there's a 
dynamically younger and /or brighter population. Also, at a 
given distance, streams with different dynamical ages (color 
scale) can be detected at different sigma levels. 



5.3 Combined Effects 

In order to combine the different effects that impact the de- 
tectability of simulated streams, explored individually in the 
previous sections, we show in Fig. QT] a plot of the total lu- 
minosity Lv of detectable satellites as a function of Dhelio, 
for all five simulated satellites at dynamical ages larger than 
4.5 Gyr and with Halo-type (bottom) and Carina-type (top) 
SFHs. This plot shows how all five satellites (S1-S5), which 
have stars observable by Gaia (and meeting our criteria 
am/zu ^ 30 per cent), can be recovered using mGC3 at 
different dynamical ages, for total luminosities in the range 
10 s -10 9 L with both SFHs and even down to 4-5 x 10 7 L Q 
for satellites including a younger stellar population (Carina- 
type SFH) for certain combinations of dynamical ages and 
orbital parameters. 

A statistical description of the efficiency of the mGC3 
method would require a more thorough exploration of the 
orbital parameter space, but in the present experiments we 
have only 5 different orbital parameter sets. Nevertheless, 
these experiments are usefull for a first rough exploration of 
the applicability of the mGC3 method in which, as shown 
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Figure 10. Height of detected excesses on subtracted pole count 
maps as a function of Dhelio- Satellites have a Carina-type SFH 
and total luminosities Lv = 2 x 10 8 L (Bottom) and Lv = 
2 X 10 9 L (Top). Detections with a > 15 are indicated with 
upward pointing arrows. The color scale represents dynamical age 
as shown by the color bar. 



in Sees. 15.11 and 15121 we find that satellites with total lumi- 
nosities as faint as ~ 5 x 10 7 L , or as dynamically old as 
~ 9 Gyr, with different orbital parameters, can be recovered. 
Furthermore, the experiments carried out illustrate a lower 
limit on the capabilities of the mGC3 method, since in our 
simulations we use only the information from trigonomet- 
ric parallaxes to obtain the distances to stars. Gaia, how- 
ever, will also provide photometric parallaxes for a much 
larger, fainter and more distant sample of stars, than those 
for which trigonometric parallaxes can be measured accu- 
rately. 



6 DISCUSSION AND CONCLUSIONS 

Motivated by the desire to unravel the fossil record of the 
formation of our G alaxy we have d e velop ed an extension of 
the GC3 method of ljohnston et all (| 19961 ) in order to apply 
this to the search for satellite remnants in the Galactic halo. 
The GC3 method is essentially a means of grouping together 
stars on the sky with the same orbital angular momentum, 
thus identifying stars along the same Great Circle. The orig- 
inal method suffers from a strong contaminating background 
in the pole count maps as demonstrated in Section [3] These 
background counts are due to the general Galactic popula- 
tion of stars and can be much reduced by adding a kinematic 
criterion to the search for stars with similar orbital angular 
momenta. This constitutes the extended mGC3 method pre- 
sented in this paper. 

We applied this new method to a mock Gaia catalogue 
containing some 3.5 x 10 8 stars that form the smooth Galac- 
tic background and stars from simulated satellite galaxies on 
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Figure 11. Total luminosity Ly of detectable satellites as a func- 
tion of the heliocentric distance, for dynamical ages greater than 
4.5 Gyr, assuming a Halo-type SFH (bottom) and a Carina-type 
SFH (top). The color scale represents dynamical age as shown by 
the color bar. 



different orbits and with two types of star formation histo- 
ries. The results show that the mGC3 method is capable 
of tracing remnants of satellites with luminosities down to 
Ly ~ 4 - 5 x 10 7 Lq for dynamical ages up to ~ 7 Gyr. 
Remnants of brighter satellites (10 s -10 9 Lq) can be recov- 
ered up to dynamical ages of ~ 10 Gyr. The method works 
well for most satellites out to a heliocentric distance of 40 
kpc. At larger distances only the brightest and/or dynami- 
cally youngest satellites can be recovered. 

Like the original GC3 method our extended version 
is limited to recovering remnants of satellites of which 
the orbits are broadly confined to a plane. This makes 
this method well suited for probing the outer halo of our 
Galaxy where the potential is more nearly spherical or ax- 
isymmetric and where dynamical timescales are long. The 
mGC3 method thus forms an important complement to the 
phase space structure characterization methods which will 
be used in the inner halo, where due to the short dynam- 
ical timescales phase space substructures are harder to de- 
tect. In these regions it will be mandatory to use very accu- 
rate measurements of integrals o f motion employing meth- 
ods such as those as propos ed in lHelmi fc de Zeeuwl (|2000| ) 
and iGomez fc Helmil (l2010h . In both cases only the use of 
the highest accuracy paralla xes (relative errors better than 
about 10 per cent, see also IGomez et all 120101 ) allows the 
recovery of phase space substructure. In contrast the mGC3 
method works with 30 per cent accurate parallax data mak- 
ing it suitable for probing larger distances. 

The Gaia mission is expected to map the immediate So- 
lar neighbourhood to very high accuracy resulting in some 
10 million stars with distances known to better than a few 
per cent. This will enable a much more accurate calibration 
of other distance indicators, in particular photometric indi- 



cators, which can then be used to extend the mGC3 method 
also to large samples beyond 40 kpc. In addition the Gaia 
mission will result in a first order smooth dynamical model 
of the Milky Way Galaxy which can be used to construct 
an accurate map of the background mGC3 pole count map. 
Subtracting a more accurate background map will further 
enhance the efficiency of the mGC3 method. 

We note here two lines of investigation that should 
be pursued to further investigate and enhance the mGC3 
method, (i) The efficiency of the mGC3 method should 
be explored more extensively by more widely sampling the 
satellite orbital parameter space. More precise limits on the 
orbital morphologies and dynamical ages, for which satel- 
lites remnants can still be recovered, can then be obtained. 
In addition, the simulation of the satellites could be made 
more sophisticated by following the suggestions made in the 
conclusions of B05. (ii) In connection with these studies it is 
interesting to investigate to what extent the requirement of 
single peaks in the pole count map for identifying remnants 
can be relaxed. This would enable us to account for the pre- 
cession of the orbital angular momentum vector and thereby 
recover a larger fraction of the stellar population of a given 
disrupted satellite. A hint of this possibility can already be 
seen in the double peaked structure of the signature of the 
satellite S2 in the map presented in Figff] 

Finally, the mGC3 method is not perfect and will cer- 
tainly pick up stars from the field population that appear 
to belong to a satellite because they are measured to have a 
compatible orbital angular momentum vector. These should 
be weeded out by making use of the fact that the different 
star formation histories for the different satellites will lead 
to different abundance patterns in their stellar populations. 
Gross distinctions (such as large overall metallicity differ- 
ences) can be made photometrically. However, eventually 
the stars identified as part of a potential satellite remnant 
should be targeted for detailed spectroscopic follow up in or- 
der to definitively 'tag' the stars to their progenitor galaxy 
and obtain accurate information on the progenitor's time of 
accretion. 

We look forward the time ahead of us when, thanks to 
the Gaia mission, complementary spectroscopic data, and 
efficient methods of characterizing substructure in phase 
space, we will enter the era of precision Galactic archae- 
ology. 
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APPENDIX A: COORDINATE 
TRANSFORMATIONS 

In the following we describe the coordinate transformation 
involved in expressing Galactocentric position r ga i and ve- 
locity v ga i in terms of the heliocentric observable quantities 
b, w, v r , Hi, Hb), latitude, longitude, parallax, radial veloc- 
ities and proper motions respectively. 

The Galactocentric position and velocity are expressed 

as 

r ga i = r + ApCf7 _1 f 

Vgai =v + v r r + uj^ 1 ^(A v fii cos fe)l + (Aj,^ 6 )bj ' ^ ' ' 

where A v = 1000 mas-pc, A v = 4.74047 yr km s _1 and 
{?, 1, b} are the unit vectors in a spherical heliocentric ref- 
erence frame. The latter are expressed in terms of the unit 
vectors in the cartesian heliocentric reference frame {x, y, z} 
as follows 

f = (cos / cos 6)x + (sin I cos b)y + (sin b)z 

I = -(sinZ)x + (cos/)y (A2) 

b = — (sinfccosZ)x — (sin b sin l)y + (cos6)z 

where x points from the Sun in the GC-Sun direction, 
y points in the direction of Galactic rotation and i points 
towards the North Galactic Pole. The relative orientation of 
the unit vectors in both reference frames is shown in Figure 

\eh 

Finally, r ga i and v ga i are 
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Figure Al. Relative orientation of the unit vectors {? , 1, b} and 
{x,y,z}, which respectively define the spherical heliocentric and 
cartesian Galactocentric reference frames used. 



r ga i = r Q + A p vj~ 1 [(cos I cosfc)x + (sinZcos&)y + (sinfr)z] 

Vgal = Vq + [v r sill 6 + AyZU -1 fib COS b]z 

+ [v r cos I cos b — A v -cu~ 1 (fii cos b sin / + fj, b sin b cos Z)]x 
+ [v r sin I cos 6 + A v vs~ (fli cos b cos I — (it sin b sin l)]y 

As can be seen from these equations, r ga i and v ga i de- 
pend upon (l,b) and the reciprocal of the parallax vj, as 
noted in Section \3. II 



